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Abstract 

Experimental cross-sections are presented for the first time for the 159 Tb(d,xn) I55 ' 157 ' 159 Dy, 155 .i 56 . lf >°Tb and 153 Gd 
nuclear reactions up to 50 MeV. The experimental data are compared with theoretical predictions of the ALICE, 
EMPIRE and TALYS nuclear reaction codes. Integral thick-target yields are also derived for the reaction products 
that have practical applications. 
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1. Introduction 

The aim of the work was manifold: 

• Very few experimental data exist in the literature 
for lanthanides, both for proton and deuteron in- 
duced reactions. 

• The quality of the database for activation cross- 
sections of deuteron induced reactions is poor, 
compared to what exists for proton induced reac- 
tion. In our laboratories a systematic study of the 
activation cross-sections of deuteron induced reac- 
tions for different applications is in progress up to 

50 Mev mm . 

• Several radionuclides of the lanthanide group are 
becoming increasingly important in the field of 
nuclear medicine (diagnostic and therapeutic ra- 
dioisotopes) SQ. 

• The prediction capability of the nuclear reaction 
model codes is still limited in case of deuteron 
induced reactions. Further improvement requires 
comparisons with experimental data. As terbium is 
monoisotopic, the targetry is simpler and cheaper 
and the cross-section results for a given activation 
product are not infected by contributions of multi- 
ple processes on the different stable target isotopes, 
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hence no need for use of enriched target technol- 
ogy- 

• Activation cross-sections on Tb are required for 
production of the medically related 159 Dy (6) and 
I57 Dy [7-9 1 and need to be compared with other 
possible production routes. 

In a literature search only two experimental cross- 
sections data sets were found, published by Due et al. 
IfTOl for production of 159 Dy, 157 Dy and I6() Tb up to 27 
MeV and by Siri et al. [ 1 1 1 for production of 157 Dy 
and 160 Tb up to 27 MeV. Mukhammedov et al. [12) 
measured integral yields of 159 Dy as a function of the 
deuteron energy up to 12 MeV for charged particle ac- 
tivation analysis. 

2. Experimental and data evaluation 

The general characteristics and procedures for irradi- 
ation, activity assessment and data evaluation (including 
estimation of uncertainties) were similar to what is dis- 
cussed in several earlier works of our group lfT3l fl4l . 
The main experimental parameters for the present study 
are summarized in Table 1. The principal methods 
used in data evaluation and the used decay data are col- 
lected in Table 2 ifBl [T31J221 and in Table 3. The good 
overlap of the re-measured excitation function for the 
27 Al(d,xn) 24 Na monitor reaction with the recommended 
values 11231 is illustrated in Fig. 1 . 
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Table 3: Investigated radionuclides, their decay properties and production routes 



Nuclide decay mode 


Half-life 


E y (keV) 


I y (%) 


Coniri billing reaction 


Q- value 










(keV) 




144.4 d 


58.0 


2.27 


liy Tb(d,2n) 


-3372.48 


e: 100 % 












157 Dj 


8.14 h 


182.424 


1.33 


15y Tb(d,4n) 


-19260.72 


e: 100% 




326.336 


93 








9.9 h 


184.564 


3.37 % 


— PTO 

Tb(d,6n) 


-35671.3 


b: 100 % 




226.918 


68.4 






iyi Dy 


6.4 h 


80.723 


11. 1 


15y Tb(d,8n) 


-51824.12 


a: 0.0094 




99.659 


10.51 






e: 99.9906 




213.754 


10.9 










254.259 


8.6 










274.673 


3.1 










389.531 


1.52 










537.225 


1.33 










1023.99 


1.09 






160 Tb 


72.3 d 


86.7877 


13.2 


15y Tb(d,p) 


4150.644 


IT: 100% 




197.0341 


5.18 










215.6452 


4.02 










298.57 83 


26.1 










879.378 


30.1 










962.311 


9.81 










966. 166 












1177.954 


14.9 






13b Tb 








'" ,y Tb(d p4n) 


-25880 05 


e: 100% 




199.19 


41 










262.54 


5.8 










296.49 


4.5 










356.38 


13.6 










422.34 


8.0 










534.29 


67 










1065.11 


10.8 










1 154.07 


10.4 










1159.03 


7.2 










1222.44 


31 










1421.67 


12.2 






155 Tb 


5.32 d 


86.55 


32.0 


^ y Tb(d,p5n) 


-32794.5 


e, 100% 




105.318 


25.1 


155 Dy decay 








148.64 


2.65 








161.29 


2.76 










163.28 


4.44 










180.08 


7.5 










340.67 


1.18 










367.36 


1.48 






l *»Gd 


240.4 d 


97.43100 


29.0 


ny Tb(d,2p6n) 


-46519.36 


e: 100 % 




103.18012 


21.1 


153 Dy decay 





When complex particles are emitted instead of individual protons and neutrons the Q-values have to be decreased by the respective binding 
energies of the compound particles: np-d, +2.2 MeV; 2np-t, +8.48 MeV; n2p- 3 He, +7.72 MeV; 2n2p-ff, +28.30 MeV 



3. Theoretical calculations 

We compared our measured cross-sections with pre- 
dictions of the ALICE-IPPE JH, EMPIRE E2 and 
TALYS 11261 nuclear reaction codes. The TALYS re- 
sults were taken from the TENDL 2011 and TENDL 
2012 nuclear reaction libraries E71 . We present both 
versions to show the improvements and the still remain- 
ing problems between the two versions of the code. In 
case of ALICE and EMPIRE we have used the codes 
modified for better description of deuteron induced re- 
actions (EMPIRE-D and ALICE-IPPE D) described in 
more detail in our previous reports Il28l[29l . 

4. Results and discussion 

4.1. Excitation functions 

The measured cross-sections for the production of the 
activation products are presented in Table 4 and Figs. 



2-8. All cross-section values of dysprosium radionu- 
clides are due to direct production via (d,xn) reactions. 
The terbium radioproducts are produced only directly 
via (d,pxn) reactions or additionally through the decay 
of the shorter-lived isobaric parent dysprosium radioiso- 
tope (cum). The investigated 153 Gd is produced di- 
rectly via (d,2pxn) reaction (including complex particle 
emission) and from the decay of parent 153 Tb radioiso- 
tope. The ground state of the produced radioisotopes 
above the direct production can be produced addition- 
ally through the internal transition of the isomeric state. 
The cross-section is marked with (m+) when the half- 
life of the isomeric state is significantly shorter com- 
paring to the half-life of the ground state and the cross- 
sections of the production of ground state were deduced 
from spectra after nearly complete decay of the isomeric 
state. 
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Table 1 : Main experimental parameters for 2 irradiation sessions 

70 




Table 2: Main parameters of data evaluation (with references) 



Figure 1 : The simultaneously measured monitor reactions for deter- 
mination of deuteron beam energy and intensity 
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4.7.7. Production of m Dy (T l/2 =144.4 d) 

According to Fig. 2 our data are systematically 
higher than the earlier results of Due et al. [ 10 1 . The 
agreement with the results of the theoretical codes is ac- 
ceptable, not counting the generally observed underes- 
timation of the high energy cross-sections by the theory. 

4.7.2. Production of l51 Dy (T 1/2 =8.14 h) 

The present and the earlier experimental data for the 
I59 Tb(d,4n) I57 Dy reaction are shown in Fig. 3. The 
agreement between all experimental data is acceptable. 
No explanation for the second bump around 50 MeV 
predicted by TENDL 2011 and 2012 was found. For 
comparison of the productions routes the cross-sections 
of Lebowitz et al. for the 159 Tb(p,3n) reaction are 
also shown. 

4.7.3. Production of l55 Dy (T l/2 =9.9 h) 

No earlier experimental data were found in the liter- 
ature. Comparison with the theory in the overlapping 
energy range shows a good agreement for EMPIRE-D, 
while the two other codes overestimate by a factor of 2 



to 5 at 50 MeV (Fig. 4). No explanations, however for 
the second bump around 90 MeV predicted by TENDL 
2011 and 2012. 

4.7.4. Production of m Tb (T 1/2 =72.3 d) 

We have only a few experimental data points at low 
energies around the maximum of the excitation function 
(Fig. 5). Our experimental data are significantly higher 
than the experimental data of Due et al [ 1 1 and a lit- 
tle lower compared to the data of Siri et al. IfTTI . The 
theoretical predictions of the two TENDL libraries are 
a factor of two low, compared to the experimental data, 
while EMPIRE-D and ALICE-D reproduce rather well 
the shape and values of our experiment. 

4.7.5. Production of l5b Tb (T l/2 =5.35 d) (m+) 

The I56 Tb is produced only directly via (d,p4n) reac- 
tion. The 156 Tb has a long-lived ground state (Ti/2=5.35 
d) and two shorter-lived isomeric states (5.3 h and 24.4 
h) decaying completely to the ground state. Our experi- 
mental data for production of ground state (Fig. 6) were 
deduced from spectra after nearly complete decay of the 
isomeric states (m+). No earlier experimental data ex- 
ist for comparison. The data in TENDL libraries show 
that the ground state is produced mostly directly. The 
contribution from the isomeric decay is small. In the 
overlapping energy range the theoretical data slightly 
overestimate the experimental results. 

4.7.6. Production of l55 Tb (T\j2—532 d) (cum) 

The measured experimental data (Fig. 7) are cumu- 
lative and are the sum of direct production and total de- 
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Figure 2: Experimental and theoretical excitation functions of the 
159 Tb(d,2n) 159 Dy reaction 

cay of the shorter lived I55 Dy (9.9 h) parent radioiso- 
tope. By referring to the excitation function of the 155 Dy 
(Fig. 4 ) the main contributor to the cumulative cross- 
section of 155 Tb is the decay of the parent 155 Dy. This 
is confirmed by the results in the TENDL libraries. The 
best result is given by the EMPIRE-D model calcula- 
tion, while ALICE-D and TENDL results overestimate 
the experimental values for the cumulative case. 

4.7.7. Production of l53 Gd (T l/2 =240.4 d) (cum) 

The experimental data for cumulative production of 
the 153 Gd are shown in Fig. 8. The experimental data 
were deduced from spectra measured after decay of the 
parent decay chain ( 153 Dy 6.4 h, 153 Tb 2.34 d). No 
earlier experimental data are available. The theoretical 
calculations display a rather large spread of the corre- 
sponding yields. For all codes the main contributions 
below 60 MeV are connected with the ff-particle emis- 
sion channels and the additional increases of the produc- 
tion cross-sections above 60 MeV are produced by the 
cumulative channels related to the decay of predeces- 
sors. TENDL-2012 underestimates distinctly the direct 
production of 153 Gd, while ALICE-D and EMPIRE-D 
strongly overestimate it. The obtained differences of the 
code results are the consequence of models used for a 
simulation of the a-particle emission. 

4.2. Integral yields 

The so called physical integral yield (instantaneous 
short irradiation) 11221 as a function of the incident 
deuteron energy was calculated from a spline fit to our 
experimental data and are shown in Fig. 9 and Fig. 
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Figure 3: Experimental and theoretical excitation functions of the 
159 Tb(d,4n) 157 Dy reaction 

10. The only measurement find in the literature is from 
Mukhamedov [12) for 160 Tb and 159 Dy. Our data are 
significantly higher than those given by Mukhamedov 
for both referred radioisotopes. 

5. Comparison of the experimental and theoretical 
results 

The comparison with the results of theoretical codes 
using global parameters shows that still large difference 
can be observed between the experimental results and 
the predictions the different codes. The new improve- 
ments for better description of the complex deuteron 
induced reactions resulted in a better description. The 
predictions for (d,pxn) reactions by of ALICE-D and 
EMPIRE-D are naturally the most successful, as the 
correction included is based more directly on the sys- 
tematic s of experimental data. The agreement with the 
2012 TENDL version is better, but the (d,p) reaction is 
still strongly under-predicted. 

6. Comparison of the production routes of 157 Dy and 

159 Dy 

Among the investigated reactions the routes leading 
to production of 157 Dy and 159 Dy have presently a prac- 
tical interest. The two radioisotopes can be produced 
via various reactions. Out of these competitive pro- 
duction routes we discuss below the low energy, light 
charged particle methods in more detail, while summa- 
rizing briefly the other routes. 
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Table 4: Measured cross-sections of the 159 Tb(d,x) 15<u57 - I55 Dy, 159 Tb(d,x) 160 - 156 ^ 155 Tb and l59 Tb(d,x) l53 Gd reactions 



E +AE (MeV) 


Cross -section a + Act (mb) 








lbu Tb 








































3.3 


1.3 


39.1 


4.5 


0.5 


0.1 






48.8 


5.5 


0.4 


0.0 










11.0 


1.2 


402.3 


45.2 


2.3 


0.3 






196.5 


22.1 


0.4 


0.1 










16.0 


1.1 


252.1 


28.4 


15.0 


1.7 






126.7 


14.2 


0.5 


0.1 










20.0 


1.0 


134.8 


15.2 


125.6 


14.1 






100.4 


11.3 


0.5 


0.1 










23.6 


0.9 


99.4 


11.2 


171.4 


19.2 






74.9 


8.4 


0.6 


0.1 










26.8 


0.8 


81.0 


9.2 


511.2 


57.4 






69.2 


7.8 


0.6 


0.1 










29.8 


0.7 


70.3 


7.9 


726.0 


81.5 






56.3 


6.3 


1.6 


0.2 










32.5 


0.7 


65.9 


7.6 


867.4 


97.4 






55.3 


6.2 


3.1 


0.4 










35.1 


0.6 


61.5 


6.9 


764.6 


85.8 






47.3 


5.3 


5.1 


0.6 










36.1 


0.6 


69.6 


7.8 


940.5 


105.6 






















37.6 


0.6 


57.0 


6.5 


654.2 


73.4 






44.4 


5.0 


9.3 


1.1 


0.4 


0.1 


0.3 


0.1 


38.2 


0.5 


63.2 


7.1 


828.0 


93.0 






















40.0 


0.5 


48.7 


5.6 


536.4 


60.2 


1.0 


0.1 


41.2 


4.6 


17.6 


2.0 


2.0 


0.2 


1.6 


0.2 


40.4 


0.5 


58.6 


6.6 


687.5 


77.2 






















42.3 


0.4 


49.2 


5.6 


433.6 


48.7 


10.9 


1.2 


37.3 


4.2 


32.2 


3.6 


12.5 


1.4 


4.3 


0.5 


42.5 


0.4 


51.7 


5.8 


558.8 


62.7 


2.2 


0.3 


















44.4 


0.4 


42.0 


4.8 


359.4 


40.3 


21.1 


2.4 


34.8 


3.9 


47.4 


5.3 


40.9 


4.5 


7.7 


0.9 


44.6 


0.4 


48.6 


5.5 


450.2 


50.5 


20.8 


2.3 


















46.5 


0.3 


41.5 


4.7 


317.0 


35.6 


70.2 


7.9 


33.2 


3.7 


65.5 


7.4 


93.6 


10.3 


11.3 


1.3 


46.6 


0.3 


44.8 


5.1 


380.8 


42.8 


86.5 


9.7 


















48.5 


0.3 


41.8 


4.8 


321.0 


36.0 


97.8 


11.0 


















48.6 


0.3 


35.8 


4.1 


260.6 


29.2 


177.0 


19.9 


31.9 


3.6 


74.5 


8.4 


195.5 


21.3 


15.4 


1.7 



Table 5: Production routes of 159 Dy and 157 Dy 



15% 


15'Dy 


Reaction 


Energy range 


Ref 


Rend inn 


Energy range 


Reference 


,W Tb(p.n) 


25-6 




8. 30 




lh9 Th(p,3n) 


30-20 


m 


l3y Tb(d,2n) 
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The I57 Dy and 159 Dy can be produced with the fol- 
lowing bombarding beams and nuclear reactions: 

• By neutron induced reactions via (n,y) on the 
neighboring stable Dy isotopes. The production 
yield is high. It requires highly enriched Dy tar- 
gets to have radionuclidic pure end-products. The 
product is of low specific activity and is not carrier 
free EH- 

• By photonuclear reactions [ 32 1 . The product has 
usually low specific activity by using (y,xn) at 
the presently used intensities and it is not carrier 
free. The specific activity could be significantly in- 
creased by using the new high intensity y-beams. 

• By high energy spallation reaction and use of elec- 
tromagnetic separators 11331 . The method assures a 
high specific activity, but the yield is low and the 
product will be expensive. It can be used in cases 
when no other method is available. Mostly in case 
of exotic and uncommon radioisotopes. 

• By low and medium energy light charged particle 



induced nuclear reactions. The production yield 
is acceptable high. In favorable cases production 
does not require highly enriched targets. The end 
product is mostly carrier free and has high specific 
activity. 

For routine production, reactor production is often the 
most favorable, if the radionuclide fulfills the require- 
ment of the applications (specific activity). Many med- 
ically interesting radionuclides can, however, be pro- 
duced only with charged particle beams. The low and 
medium energy, high power (beams up to 1 mA) cy- 
clotrons are also becoming more competitive for pro- 
duction of those radioisotopes, which presently pro- 
duced in nuclear reactors. Photonuclear and spalla- 
tion based methods are still not common, these methods 
need improvements on the beam intensity. The main 
production routes for 157 Dy and 159 Dy via light charged 
particle beam are summarized in Table 5 Q [8] [10] [TT] 
l30l . The excitation functions of the reactions are shown 
in Fig. 1 1 and 12. 

Taking into account that for many of the involved ex- 
citation function no experimental data exist, the excita- 
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Figure 4: Experimental and theoretical excitation functions of the 
159 Tb(d,6n) 155 Dy reaction 



Figure 5: Experimental and theoretical excitation functions of the 
159 Tb(d,p) 160 Tb reaction 



tion functions are based on the adjusted theoretical pre- 
dictions of the TENDL 2012 library. The differences 
in the different routes are so high that conclusion is not 
affected by the known poor prediction capability of the 
theoretical results. 

6.1. Production routes of 159 Dy 

According to Table 5 and Fig. 1 1 the most favorable 
route for production of 159 Dy at a low energy cyclotron 
(k = 20) is the 159 Tb(p,n) process and at higher energies 
(k = 35) the 159 Tb(d,2n) reaction. The yield of the (d,2n) 
reaction is, like in many cases we discussed earlier, sig- 
nificantly higher than the yield of (p,n). According to 
Fig. 1 1 the cross-sections of the 3 He induced reactions 
are very small, similar to cross-sections of the (a,n) re- 
action. The (a,2n) process requires a 30 MeV cyclotron 
and highly enriched target. The cross-section is high but 
has lower yield comparing to the (d,2n) reaction due to 
the shorter range of the alpha particle. The 159 Dy can 
also be produced with high radionuclide purity on na 'Gd 
by both ( 3 He,xn) and alpha particle induced reactions. 
In both cases the yield is low. 

6.2. Production routes of 157 Dy 

According to the Fig. 12 at a low energy cyclotron 
only the low yield 154 Gd(a,n) reaction is available. At 
commercial 30-35 MeV FT-cyclotrons the 159 Tb(p,3n) 
reaction results in high production yields. In the opti- 
mum (30-20 MeV) energy range the amount of the si- 
multaneously produced 159 Dy is small due to the low 
cross-sections of the (p,n) in this higher energy range 



(see Fig. 12). The energy range for the (d,4n) reac- 
tion should be 50-30 MeV to minimize the amount of 
159 Dy produced simultaneously via (d,2n). Accelera- 
tors capable of practical implementation of the (d,4n) 
reaction are however very rare. The 155 Gd(a,2n) re- 
quires highly enriched targets and cyclotrons with 30 
MeV alpha beams, and the yield is significantly lower 
compared to the 159 Tb(p,3n) reaction. To get high ra- 
dionuclide purity highly enriched 155 Gd and 156 Gd tar- 
gets should be used in case of 3 He-particle induced reac- 
tions. The yield however is very low and the 3He beam 
is expensive even having recovery system for the gas in 
the ion source. 



7. Summary and conclusion 

The present work provides new data for the excita- 
tion function of deuteron induced nuclear reactions on 
natural terbium, which is actually 159 Tb. The experi- 
mental results for the produced radioisotopes were com- 
pared with the former literature values as well as with 
the results of theoretical model calculations by using 
the codes ALICE-D, EMPIRE-D and TALYS 1.4 (the 
data taken from the TENDL 2011 and 2012 libraries). 
Our new measurements both verified or improved the 
former results and gave an extension to them. Compar- 
ison with the results of model codes could validate the 
recent improvements of those, as well as provide new 
basis for further improvements. As practical results, 
physical yield curves were also calculated for medical 
and other practival use. Serving for the present demands 
of radioisotopes research a detailed comparison is given 
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Figure 6: Experimental and theoretical excitation functions of the 
159 Tb(d,p4n) 156 «Tb(m+) reaction 
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for the possible routes for production of 2 important ra- 
dioisotopes 157 Dy and 159 Dy. 



Figure 7: Experimental and theoretical excitation functions of the 
159 Tb(d,x) 155 Tb(cum) reaction 
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Figure 8: Experimental and theoretical excitation functions of the 
159 Tb(d,x) 153 Gd(cum) reaction 
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Figure 9: Thick target yields for radionuclides of dysprosium pro- 
duced by deuteron irradiation on I59 Tb and comparison with the avail- 
able literature values 



Figure 11: Comparison of excitation functions for production of 
159 Dy 
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Figure 10: Thick target yields for radionuclides of terbium and 
gadolinium produced by deuteron irradiation on 159 Tb and compar- 
ison with the available literature values 



Figure 12: Comparison of excitation functions for production of 
157 Dy 
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